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A commentary on 

Commentary on "Force-field functor the- 
ory: classical force-fields which reproduce 
equilibrium quantum distributions." 

Vaia, R. (2013). Front. Chem.l:34. doi: 10. 
3389/fchem.2013.00034 

Dr. Ruggero Vaia has made several claims 
(Vaia, 2013) about our recent submission 
entitled Force-field functor theory: classi- 
cal force-fields which reproduce equilibrium 
quantum distributions (Babbush et al, 
2013) which we would like to address. Dr. 
Vaia points out that the idea of an effec- 
tive potential which reproduces quantum 
distributions is a well-known observation. 
We agree and point out that our paper 
contains many citations to the relevant lit- 
erature. The contribution we make is to 
describe a new perspective on how the 
effective potential might be obtained in a 
spirit similar to Density Functional Theory 
(DFT). To make progress on a DFT- 
like theory, a proof of the existence and 
uniqueness of force-field functors is nec- 
essary. The main goal of our manuscript is 
to develop such proofs and to demonstrate 
and test the idea of functor on various test 
problems. 

After establishing uniqueness and exis- 
tence proofs, we introduce an intention- 
ally simple example of a functor in order 
to illustrate the utility of our broader 
approach. In particular, we use the Jensen- 
Peierls inequality to construct a variational 
approximation in which the mapping is 
linear. This procedure closely parallels the 
effective potential of Feynman and Hibbs 
(1965) in which the same approximation 



is used to model the effective potential as a 
linear convolution of the physical potential 

with a Gaussian of variance = ^ m. 
Dr. Vaia correctly shows that for a one- 
dimensional system, our analytical approx- 
imation reduces to a Gaussian smearing 
of the potential with exactly twice the 
width of the Feynman-Hibbs approxima- 
tion. This is because the Feynman-Hibbs 
potential is a centroid effective poten- 
tial which does not attempt to reproduce 
the correct particle distribution. While we 
make reference to this in our paper (we 
wrote "our [linear mapping] can be imag- 
ined as a Gaussian smearing of y(q)"), 
we do not discuss the explicit form of 
the mapping in one-dimension as this 
obscures the more general point we are 
trying to make. 

Linearity is not an essential feature of 
our approach; rather, it is an approxima- 
tion which underlies a particular functor 
that we construct as a pedagogical exam- 
ple of how the functor approach could be 
used. In fact, the procedure we demon- 
strate to invert the linear mapping is sig- 
nificantly more accurate than the explicit 
Gaussian convolution which Dr. Vaia dis- 
cusses. We used the bijective nature of 
the mapping (implied by our uniqueness 
proof) to empirically fit the optimal linear 
functor. This was accomplished by com- 
puting the exact effective potentials for 
1000 random physical potentials and then 
using a least squares procedure with cross- 
validation to determine optimal matrix 
elements which map basis vectors of the 
physical potential to basis vectors of the 
effective potential. This cross-validation 
approach allowed us to test the accuracy 



of our method. We show that our functor 
works surprising well (especially consid- 
ering its simplicity) and demonstrate its 
impressive performance on independent 
test problems including the simulation of 
liquid hydrogen at temperatures as low as 
14 K. Though the error in this approxima- 
tion grows as temperature decreases, the 
same is true of all path integral approaches 
since the number of required integration 
time slices grows without bound as tem- 
perature approaches zero. Degrading in 
quality at T = 0 does not render path inte- 
gral molecular dynamics, or our technique 
useless. 

Dr. Vaia goes on to compare the linear 
approximation to the Feynman-Kleinert 
and Giachetti-Tognetti effective potentials 
and mentions that those approximations 
are more robust. While this is true, we are 
not trying to demonstrate the superiority 
of a particular explicit form for the func- 
tor under any particular approximation and 
we are certainly not interested in comparing 
the performance of this simple example to 
other paradigm methods in one dimension. 
The linear functor is for our theory analo- 
gous to what Local Density Approximation 
(LDA) is for DFT: a simplistic but illustra- 
tive example of the type of theory which 
could emerge from our uniqueness the- 
orem. We have provided ample evidence 
for a new approach to obtaining effective 
potentials, in a similar spirit to DFT. 
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